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Abstract 

A program for computing pure angular momentum coefficients in relativistic atomic 
structure for any scalar one- and two-particle operator is presented. The program, writ- 
ten in Fortran 90/95 and based on techniques of second quantization, irreducible tensorial 
operators, quasispin and the theory of angular momentum, is intended to replace existing 
angular coefficient modules from GRASP92. The new module uses a different decom- 
position of the coefficients as sums of products of pure angular momentum coefficients, 
which depend only on the tensor rank of the interaction but not on its details, with ef- 
fective interaction strengths of specific interactions. This saves memory and reduces the 
computational cost of big calculations signficantly. 



PROGRAM SUMMARY 



Title of program: ANCO 
Catalogue number: ADOO 

Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland. 
Users may obtain the program also by down-loading the tar-file ratip-anco . tar from our 
home page at the University of Kassel 



Program Summary URL: http://cpc.cs.qub.ac.uk/summaries/ADOO 
Licensing provisions: None. 

Computer for which the program is designed and has been tested: 

IBM RS 6000, PC Pentium II. 

Installations: University of Kassel (Germany). 

Operating systems: IBM AIX 4.1.2+, Linux 6.1+. 

Program language used in the new version: ANSI standard Fortran 90/95. 
Memory required to execute with typical data: 100 kB. 

No. of bits in a word: All real variables are parametrized by a selected kind parameter. 
Currently this is set to double precision for consistency with other components of the RATIP 
package [1]. 

Distribution format: Compressed tar file. On a UNIX (or compatible) workstation, the 
command uncompress uncompress this file and the command tar -xvf ratip_anco .tar 
reconstructs the file structure. 

Keywords: atomic many-body perturbation theory, complex atom, configuration interaction, 
effective Hamiltonian, energy level, Racah algebra, reduced coefficients of fractional parent- 
age, reduced matrix element, relativistic, second quantization, standard unit tensors, tensor 
operators, 9/2-subshell. 

Nature of the physical problem: 

The matrix elements of a one-electron tensor operator A k of rank k with respect to a set 
of configuration state functions \~fiJi) can be written J2 a b^ij( a ^) (a\A k \b) where the angular 
coefficients t k j(ab) are independent of the operator A k , i,j are CSF labels and a, b run over 
the relevant interacting orbital labels. Similarly, the matrix elements of the Dirac-Coulomb 
Hamiltonian can be written in the form Y^ab^iji^) { a \HD\b) + YkJ2 abed v k j(abcd) X k (abed) , 
where Hp is the one-electron Dirac Hamiltonian operator, with tensor rank zero, v k j(abcd) 




http:/ / www.physik.uni-kassel.de /fritzsche / programs.html ) . 
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are pure angular momentum coefficients for two-electron interactions, and X (abed) denotes 
an effective interaction strength for the two electron interaction. The effective interaction 
strengths for Coulomb and Breit interaction have different selection rules and make use of 
subsets of the full set of coefficients vfj (abed) . 

Such matrix elements are required for the theoretical determination of atomic energy levels, 
orbitals and radiative transition data in relativistic atomic structure theory. The code is 
intended for use in configuration interaction or multiconfiguration Dirac-Fock calculations [2], 
or for calculation of matrix elements of the effective Hamiltonian in many-body perturbation 
theory [3]. 

Method of solution: 

A combination of second quantization and quasispin methods with the theory of angular 
momentum and irreducible tensor operators leads to a more efficient evaluation of (many- 
particle) matrix elements and to faster execution of the code [4] . 

Restrictions on the complexity of the problem: 

Tables of reduced matrix elements of the tensor operators a^ q ^ and W^ kqk ^ are provided for 
(nj) with j = 1/2,3/2,5/2,7/2, and 9/2. Users wishing to extend the tables must provide 
the necessary data. 

Typical running time: 

3.5 seconds for both examples on a 450 MHz Pentum III processor. 
Unusual features of the program: 

The program is designed for large-scale atomic structure calculations and its computational 
cost is less than that of the corresponding angular modules of GRASP92. The present 
version of the program generates pure angular momentum coefficients i° (ab) and vfj (abed) , 
but coefficients t^(ab) with k > have not been enabled. An option is provided for generating 
coefficients compatible with existing GRASP92. 

Configurational states with any distribution of electrons in shells with j < 9/2 are allowed. 
This pemit user to take into account the single, double, triple excitations form open d— and 
/- shells for the systematic MCDF studies of heavy and superheavy elementds (Z > 95). 
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LONG WRITE-UP 



1 Introduction 

The improved accuracy of modern experiments challenges theorists to match or exceed ex- 
perimental precision. Models of many-electron atoms and ions require both relativistic and 
correlation effects to be taken into account; this can be done, for example, by using various 
versions of perturbation theory, the configuration interaction method, the multiconfiguration 
Hartree-Fock method d or the multiconfiguration Dirac-Fock method [2]. 

The evaluation of matrices of one- and two-electron operators for many-electron states in 
jj-coupling is customarily done by expressing each matrix element as a sum of products 
of angular coefficients and radial integrals. This strategy, based on earlier work by Fano, 
was adopted for the MCP program 11] for evaluating angular coefficients for the Coulomb 
interaction and the related MCBP program [5] for evaluating angular coefficients for the Breit 
interaction. Whilst this strategy was adequate on the relatively low-powered computers of 
the 1970s, it is now possible and desirable to use very large configuration sets which require 
a more efficient strategy. The new program attains this objective by building up the angular 
coefficients for one- and two-electron interactions from a relatively small number of common 
spin-angular parts in the manner of (HIIZIIHI- 

The theoretical background is presented in Section 2, program organization in section 3, 
testing and timing studies in section 4 and simple test problems in Section 5. 

2 Theoretical background 
2.1 Dirac— Coulomb Hamiltonian 

As usual, multiconfiguration self-consistent-field calculations in relativistic atomic theory are 
based on the Dirac-Coulomb hamiltonian: 



and the second term in Q represents the Coulomb interaction of pairs of electrons. In 
equation (J2J) c denotes the speed of light; oc(i) and f3{i) are 4x4 Dirac matrices for the 
ith electron; r{ and Pj are the radial coordinate of the ith electron and its (3-) momentum, 
respectively. The first two terms of (J2J comprise the Dirac kinetic energy operator. 

Both configuration interaction and multiconfiguration Dirac-Hartree-Fock calculations re- 
quire the matrix elements of Hoc with respect to a basis of n-electron configurational states 





Ho{i) is the Dirac one-particle operator, 




(2) 
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labelled |7 r </r), ^ = 1,2,3,... The present program computes the pure angular coefficients 
t k s (ab) (only for k = in the present version) and v k s (abcd)) in the expression 

{ lr Jr\H DC \lsJs) = E {t°rs(ab) (a\H D \b)+^J24s\cibcd)X^(abcd)^ > ( 3 ) 



u6 



k cd 



where X^ k \abcd) is the effective interaction strength of the two-electron interaction with 
respect to the orbitals concerned. Here in (jHJ), the interaction is the Coulomb potential only, 
where 



X {k \abcd) 



(4) 



(-in 

\\C^\\n c l c j c ) {n b l b j b \\C^\\n d l d j d )R k (n a l a j a n b l b j b n c l c j c n d l d j d ) 



but the formalism can also be used for the Breit interaction, with a different definition of 
X^ (abed) and different selection rules. Equation Q is a rearrangement of the formula used 
in GRASP92 9, eqn. (3.10)], 



(lrMH DC \lsJs) = £ T rs (ab)I(ab) + ]T ]T V^\abcd) R^(abcd) 



(5) 



ab 



k cd 



where 



and 



I(ab) = [a\H D \b 



R^(abcd) = (a&|r< fc /r> fc+1 |cd) 

is a relativistic Slater integral over orbitals a,b,c,d (in the usual notation). Whilst the 
one-electron part is the same, equation © expands the effective interaction strength as a 
traditional sum of Slater integrals, and Vrs\abcd) therefore differs from v^\abcd). The latter 
can be used for any two-electron interaction, and it is no longer necessary to treat Coulomb 
and Breit interactions on a different footing in the manner of GRASP92. Other infractions 
such as the lowest-order normal mass shift and the specific mass shift eqns. (3.12-3.14)] 
can be handled within the same scheme. The way in which the coefficients are built up is 
described below. 



2.2 One— particle operators 

The matrix elements of a one-particle scalar operator between configuration state func- 
tions with u open shells can be expressed as a sum over one-electron contributions 



where 

(^(J)\\F(n i Ki 1 n jKj )\\^(J' 

= (-l) A+1 v ^j7TTi?(j l , Jj ,A frra ,A fcet ) 5(Ki, Kj ) ( 



II i It i. 



f(0) 



(6) 
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x <! S(m,nj) [jf 1 aiQiJj 



Jg k) v _(? 3i) 

"l/2 X a -l/2 



(0) 



+ 



(1 - S(rii,nj)) h'f* aiQjJi 



,(<?ii) 

l l/2 



if' "iQi-h) (if 3 "A ; r 7 / 



l -l/2 



All states are defined in jj-coupling. (^„ m (J) || and HV'u 6 * ("^0) are respectively bra and 
ket functions with u open subshells, k = (2j + 1)(Z — j), ^n-jKj /(°) njKj) is the one- 
electron reduced matrix element of the operator , A bra = (Jj, Jj, , Jj') bra and A fce * = 
(Jj, Jj, Jj/, Jj') ket denote the respective sets of active subshell angular momenta. The opera- 
tors a^, are second quantization operators in quasispin space of rank q = 1/2. The operator 



Mi) 

6 l/2 mj 



os creates electrons with angular momentum quantum numbers j, rrij and its 



conjugate a 



(qj) 

■l/2mj 



;(J) 



(— 1) J mj a^l^. annihilates electrons with the same quantum 



numbers j, rrij in a given subshell. 
2.2.1 Recoupling matrix 

The recoupling matrix R (ji , jj , A frra , A ket ^j in © is particularly simple. It is either a product 
of delta functions [BJ eqn. (18)] when combination of delta functions and 

6j— coefficients eqn. (22)] when 7^ 



2.2.2 Matrix elements of irreducible tensor operators 

By applying the Wigner-Eckart theorem in quasispin space we obtain the submatrix elements 
of operators of type am, in the form ^U] 



(j N aQJ 



Mi) 



~[Q\ 



-1/2 



j N a'Q'J' 



Of 1/2 Q 

M' Q m q M Q 



Mi) 



j a'Q'J'), 



(7) 



where we have used the conventional shorthand notation (2k + 1) • ... = [A;, ...], and the last 
factor is a reduced coefficient of fractional parentage. The submatrix elements of the simplest 
compound tensor operator of type 

nj N aQJ 

= E iQr 1/2 



m q i m q 2 



(k 3 ) 



uses 

■N' 1 r^i ji 
nj aQJ 



m ql m q2 m q 

-11 — y 

x (nj aQJ hrfa^l nj a'Q'J') . 



Q fcq Q 

M' Q m q M Q 



(8) 



where [nj aQJ 



nj a'Q'J') denotes the reduced matrix element of the tensor oper- 



ator W^ kqk ^ (nj,nj) = a^ q ^ x a^ q ^ q J ^ in quasi-spin space. In terms of the fully reduced 
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coefficients of fractional parentage (j aQJ 

nj a'Q'J' 



j a'Q'J' ) , we find 



(_l)«+J+«'+J-+*,+*i[^ |fci ]V2 ^ 



q q k q 

a"Q"J" I Q Q Q 



J' J J" 



• H n>i t»\ I ■ a nil j s 
ja Q J )[ja Q J 



Mi) 



j a'Q'J 1 



(9) 



This construction has the advantage that the completely reduced matrix elements on the right 
hand side of (J7J) and (jHJ) are independent of the occupation number of the shell, and so reduces 
requirements of storage in comparison with earlier work. These formulae are evaluated in the 
module rabsjrcfp [TT] . 

The phase factor arises from the reordering needed to match the recoupled creation and 
annihilation operators in the bra and ket vectors. We have 



when ; otherwise 



0. 



b-l 



(10) 



(11) 



where N r is the occupation number of subshell r, a = min{i, j}, and b = max{i, j}. 

2.2.3 The one electron submatrix element 

It only remains to define the one-electron interaction matrix element 

in The only operator required in this implementation is the matrix element of the Dirac 
operator, a tensor operator of rank zero, 



. TbiKi 



(12) 



where I(riiliji,riiliji) is defined by |12l eqn. (22)]. The Dirac kinetic energy operator, denoted 
by T in jUJ eqn. (3.13)], can be obtained from this by setting the nuclear charge Z = 0. The 
coefficients T rs (ab) in |2J) can now be identified by inserting 1)12(1 in @. Meanwhile the pure 
angular coefficients t®. s (ab) for one-electron operators can be identified by inserting 



(n iKi \\fM 



m 
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2.3 Two-particle operators 

According to jS], the matrix element of any two-particle scalar operator G^ kk0 ^ between 
configuration state functions with u open shells, can be written 



E 



(^u a (J) GiniK^njKj^i'Ki^rij'Kf) tp ket (J')) (13) 



where 



(ifj b u ra (J) \\G (ni^njKj^i'Ki^nfKf) \\ip ket (J')) 

fcl2 

xT (n i ji,n j j j ,n il j i >,n jl j jl ,K hra , A ket , 3, r) (ii , jj , ji' , j/ , A 6 ™ , A fce * , r) , 

r specifies the recoupling scheme required for each matrix element and 5, when required, 
specifies the coupling scheme of the tensor operators defining each matrix element. The 
operator G^ kk0 ^ couples tensor operators of rank k for each electron to give an overall scalar 
operator. 

From (|13[) we see that the matrix element of any two-particle operator can be written as 
a sum over all possible sets of active shell quantum numbers The 
systematic analysis of jH] aims to minimize the computation needed in this expansion. The 
parameter distributions are presented in Table ^ Note that for distributions 2-5 and 19- 
42 the subshell labels are ordered so that a < f3 < 7 < 5, while for distributions 6-18 no 
conditions upon the ordering are imposed. We discuss these structures in more detail below. 



2.3.1 Recoupling matrix 

The recoupling coefficients defined in 012121 did not reduce the recoupling coefficients to their 
simplest forms but relied on the analysis module of the NJSYM package (later NJGRAF) 
to perform the reduction mechanically. The analysis jHj leads to simpler forms denoted by 
R (ji , jj , , jji , A bra , A ket ,Fj. In the case of one interacting shell R (ji , jj , jV , jji , A bra , A ket , T^j 
reduces to delta functions [HI eqn. (18)]. For two, three and four interacting shells, the recou- 
pling coefficients are given by 6, eqns. (22), (23) and (24)], replacing I, L by j, J respectively. 
The recoupling parameters T for each distribution can be found in Tabled 



2.3.2 Matrix elements of irreducible tensor operators 



The expressions T yiiji,njjj,riiijii ,nj>jj' 1 A bra , A ket , H,rj are matrix elements of standard 
subshell creation/annihilation operators 



W 



,(») 



a rn q \ X a m q \ 



(14) 
(15) 
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Table 1: Scheme of the definitions for matrix elements of any two-particle operator. The 
operators a, W, aW and Wa defined in ((Til) - (fTS)) act on the indicated subshells. 



No. 


i 


3 


i' 


f 


r 




a 


P 


7 


5 


<P 


1. 


a 


a 


a 


a 


- 


3a, 


k 


WW 


- 


- 


- 


— 


















- 


ja 




W 


- 


- 


— 


— 








2. 


a 


P 


a 


P 


k 


3a, 


30 


W 


W 


- 


— 











3. 


P 


a 


P 


a 


k 


ja, 


30 


W 


W 


- 


- 











4. 


a 


P 


P 


a 


ki2 


ja, 


30 


W 


w 


- 


- 











5. 


P 


a 


a 


P 


k\2 


3a, 


30 


W 


w 


- 


- 











6. 


a 


a 


P 


P 


k\2 


3a, 


30 


W 


w 


- 


- 


ja 


+ 30 


+ k 12 




7. 


P 


a 


a 


a 


.?> 


3a, 


ki2 


aW 


a 


- 


- 


k- 


^12 






8. 


a 


P 


a 


a 


jp 


3a, 


k\2 


aW 


a 


- 


- 


k 








9. 


p 


P 


P 


a 


ja 


30, 


k\2 


a 


Wa 


- 


- 


k- 


^12 






10. 


p 


P 


a 


P 


3a 


30, 


k\2 


a 


Wa 


- 


- 


k 








11. 


P 


7 


a 


7 


3a, jp, k 






a 


a 


W 


- 


1 - 


h ja 4 


-30 ~ 


k 


12. 


7 


P 


7 


a 


3a, 313, k 


- 




a 


a 


W 


- 


1 - 


h ja 4 


-30 ~ 


k 


13. 


7 


P 


Q 


7 


3a,jf3,kl2 


- 




a 


a 


w 


- 


1 - 


h ja 4 


-30 - 


k\2 


14. 


P 


7 


7 


a 


3a,jf3, kl2 


- 




a 


a 


w 


- 


1 - 


h ja 4 


-30 - 


k\2 


15. 


7 


7 


Q 


P 


ja, 30, k\2 


— 




a 


a 


w 


- 


ja 


+ 3i 


+ kn 




16. 


7 


7 


P 


a 


3a, 30, k\2 


— 




a 


a 


w 


- 


ja 


+ ji 






17. 


a 




7 


7 


3 a, 30, k\2 


— 




a 


a 


w 


- 


ja 


+ 30 


+ ki2 




18. 





Q 


7 


7 


3a, 30, k\2 


- 




a 


a 


w 


- 


30 


+ ji 






19. 


a 


P 


7 


S 


3a, 3 0,3 1,35 


— 




a 


a 


a. 


a 


ja 


+ 30 


+ k 12 




20. 


P 


a 


7 


S 


3a, 3 0,3 1,35 


- 




a 


a 


a 


a 


30 


+ ji 


+ k 12 




21. 


a 


P 


S 


7 


ja, 30, 3 1,35 


- 




a 


a 


a 


a 


30 + ji 






22. 


P 


a 


S 


7 


ja,j0,ji,js 


- 




a 


a 


a. 


a 


30 + h 






23. 


7 


5 


a 


P 


ja,j0,ji,j5 


- 




a 


a 


a 


a 


ja 


+ 36 


f k 12 




24. 


7 


5 


P 


a 


ja,j0,ji,js 


- 




a 


a 


a 


a 


ja 


+ js 


f h2 




25. 


5 


7 


a 


P 


ja,j0,ji,j5 


— 




a 


a 


a. 


a 


ja 


+ js 






26. 


5 


7 


P 


a 


ja,j0,ji,js 


- 




a 


a 


a 


a 


ja 


+ js 






27. 


a 


7 


P 


5 


ja,j0,ji,js 


- 




a 


a 


a 


a 











28. 


a 


7 


S 


P 


ja,j0,ji,j5 


- 




a 


a 


a. 


a 











29. 


7 


a 


S 


P 


ja,j0,ji,js 


— 




a 


a 


a 


a 











30. 


7 


a 


P 


S 


ja,j0,ji,js 


- 




a 


a 


a 


a 











31. 


/? 


6 


a 


1 


ja,j0,ji,js 


- 




a 


a 


a. 


a 


ja 


+ 30 


+ ji - 


V36 


32. 


8 


P 


7 


a 


ja,j0,ji,js 


— 




a 


a 


a 


a 


ja 


+ 30 


+ 3i - 


V36 


33. 


P 


5 


7 


a 


ja,j0,ji,j5 






a 


a 


a 


a 


ja 


+ 30 


+ ji - 


V36 


34. 


S 


P 


a 


1 


ja, jp, 3i, 36 






a 


a 


a. 


a 


ja 


+ 30 


+ 3i- 


V3& 


35. 


a 


S 


P 


7 


ja,j0,ji,js 






a 


a 


a 


a 


1 - 


Hj 7 4 


36 ~ 


k 


36. 


S 


a 


7 


P 


ja,j0,ji,js 






a 


a 


a 


a 


1 - 


hj 7 4 


36 - 


k 


37. 


a 


S 


7 


P 


ja,j0,ji,js 






a 


a 


a 


a 


1 - 


hj 7 4 


36 - 


ku 


38. 


S 


a 




7 


ja,j0,ji,js 






a 


a 


a 


a 


1 - 


Hj 7 4 


js - 


k\2 


39. 


P 


7 


n 


5 


ja,j0,ji,js 






a 


a 


a 


a 


1 - 


h ja 4 


-30 - 


k 


40. 


7 


P 




a 


ja,j0,ji,js 






a 


a 


a 


a 


1 - 


h ja 4 


-30 ~ 


k 


41. 


P 


7 


6 


a 


ja,j0,j 7 ,js 






a 


a 


a. 


a 


1 - 


h ja 4 


-30 - 


k\2 


42. 


7 


P 


a 


S 


ja,j0,ji,jg 






a 


a 


a 


a 


1 - 


h ja 4 


-30 ~ 


ki2 
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aW 
Wa 
WW 



x 



,(»■) 



n to) x n (<3i) 



(fcia)' 



m 9 2 



x a. 



™ 9 3 



(*a) 



(A;) 



"m,3 U "iq4 



(0) 



(16) 
(17) 
(18) 



The creation and annihilation operators in (|14|) - (|18|) refer to a single subshell. The evaluation 
of the submatrix elements of operators of type a (|14[) and the simplest compound tensor 
operator of type W (fTo]) was explained in section 2.1.2. For types lfH))l -lfT5 |) . we use the 
formula 



nj N aQJ 



U {kl) (nj) x V^(nj) 



(fc) 



■N' i rM t/ 

nj aQJ 



(19) 



{ _ 1)J+ j /+k[k]1/2 ^ 

a"Q"J" 



k\ hi k 
J' J J" 



x (nj N aQJ \\u^(nj)\\ nj N " a"Q"J") (nj N " a"Q"J" \v^(nj)\ nj N ' a'Q'J' 

where U^ kl ^ (nj) , y( fc2 ) (nj) are either of type (|14(1 or type (j!5j) . The occupation number 
N" is defined by second quantization operators occurring in (nj) and V (k ^ (nj). The 
module rabs_rcfp |llj performs the evaluation of these formulae. 

2.3.3 Phase factors A 

These arise from the reordering necessary to match the recoupled creation and annihilation 
operators in bra and ket vectors contributing to the matrix element. For each of the cases 
considered in Table ^ we find 



Cases 1-6: 



Cases 7-18: 



0. 



3-1 

A = l + ^iV r , 



(20) 



(21) 



where N r is the occupation number of subshell r. If a < (3, then i = a, j = (3, and if a > /3, 
then i = (3, j = a. 

Cases 19-42: 



P-l (5-1 

A = E N k + E N k- 

k=a k=-y 

2.3.4 The coefficients B' (niliji,njljjj,ni/li'jii,njdj'jj',S) 

The effective interaction strength of order A; of a two-electron operator, is denoted by 



(22) 



Jkk) 



ni'li'jiinjdj'jj' 
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in [Jj and by 



X k (riikji , rijljjj , riidi'ji' , nj'lfjji 



in ^2]. The coefficients 0' (riikji,rijljjj,Tiidirji',rijdjrjjr,B,) for the different cases tabulated 
in Table^) have different multiplicative factors defined as follows: 

Case 1 : Single subshell (aaaa) 

®JJa {nalajat n a l a j a , n a l a j a , n a l a j a , 



\k\ I (n a l a j a Tl a l a j a 



Jkk) 



(23) 



and 



®IIb {nalajai n a l a j a , n a lajai n a lajai ') 



( 1) [jet] ^ (TL ct l ct j a Tl a l a j, 



(kk) 



Lja) 8 (kn,0) 



(24) 



Cases 2, 3, 11, 12, 27, 29, 31, 32, 35, 36, 39, 40: Subshell assignments a/3a/3, papa, p-yaj, 
7/?7a, «7/3<5, ~jab~P, /3<5cry, 5P^a, aSP^f, 5ajP, P^/a5, ^P5a 



©' (nikjijnjljjjjrii'li'ji'jnj'lj'jj'jE) 



{~ l Y \ [k] 1/2 (nikjiTijljjj # (fcfc) riidi'ji'nflj'jj^ 5 (k 12 , k) . 



(25) 



Cases 6, 15-26: Subshell assignments aaPP, 77a/?, 77/3a, a/?77, /3a77, a/?7<5, Padj, aPdj, 
Pwy5, ^/5aP, S^fPa, 'ySPa, 5^aP 



0' (niliji,njljjj,ni'li'ji',nj'lj'jj',E) 



1 r, ,1/2 J ji ji' k 



(-l) 1+fc+ %[^ , . . , 

2 { 3-y h k\2 



(nikjiTijljjj 5 (fcfe) ni'k/ji/nj'lj'jjf) . (26) 



Cases 4, 5, 7, 8, 9, 10, 13, 14, 28, 30, 33, 34, 37, 38, 41, 42: Subshell arrangements aPPa, 
PaaP, Paaa, aPaa, PPPa, PPaP, jPa^f, P^a, a^bP, jaP5, P8^a, 5/?«7, a5^P, 5aPj, 
P~/5a, ^Pa5 



0' {nikji,njljjj,n i di>ji>,njiljijji,'E) 

= {- l Y\ [k\2\ l/2 \ .* Ji \ (nikjiTijljjj 

2 I Jj 3j> Ki9 I v 



{kk) 



riidiijiirijdjijji 



(27) 



The phase factors ip in expressions (|25|) - (|27|) are defined in column </j of Table ^ This 
construction exploits the common tensorial structure of any scalar two-electron operators 
as the Coulomb, Breit and Gaunt interactions ^2]) and exploits this similarity to simplify 
the calculation of spin-angular coefficients. The relativistic jj-coupling expressions for the 
effective interaction strength of the Coulomb interaction j!21 eqn. (86)] is 



(riikjinjl 



Jkk) 



nidiiji/Tljdjijji 



(28) 



{-\) k (nikji\C^\niditjit){njljjj\C^ 
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We can now identify the coefficients V*. {abed) of © by substituting the results above in . 
The same construction can be used for the Gaunt interaction (the leading part of the magnetic 
Breit interaction) |121 ean. (91)] and for the full transverse Breit interaction |12l ean. (101)], 
although the selection rules and the effective interaction strengths corresponding to (|27|) are, 
of course, different. The pure angular coefficients v^ s (abcd) for two-electron operators are 
the same for all these operators since the v^ s (abcd) can be identified by inserting 



(nihjirijljjj £r (fefc) nidi'jitrijdj'jj 



in <H3J). 



3 Program organization 

3.1 Overview of the program 

The program ANCO constructs the pure angular coefficients t® s (ab) for one-electron operators 
and the v^ s (abcd) coefficients contributing to matrix elements of the Dirac-Coulomb-Breit 
Hamiltonian. The coefficients T rs (ab) and V^ s (abcd) used in GRASP92 and earlier version of 
the system are available as an option. The new format generates what we have called "pure" 
angular momentum coefficients which can be used unchanged with any one-particle tensor 
operator of rank 0, and any two-particle interaction. The Coulomb and Breit interactions 
use different subsets of the complete set of v^ s {abcd) coefficients, which are selected automat- 
ically when multiplying by the relevant effective interaction strengths to complete the matrix 
element calculation. The MCP and MCBP modules of GRASP92 calculated the full matrix 
elements for each of these subsets, so that the new formulation reduces the computational 
overheads and the memory requirements, which renders ANCO more suitable for large scale 
problems. 

There are two new modules rabs_recoupling and rabs_anco for extracting spin-angular 
coefficients relating to formula (|13j) . The module rabs_recoupling evaluates recoupling 
coefficients R (j u j j} A bra , A fcet ) and R(j i ,j j ,j i/ ,j f ,A bra ,A ket ,r\ as described in :01, mod- 
ule rabs_rcfp evaluates the T (niji,njjj,riiijir,nj>jji, A fera ,A fcei ,H,r) whilst rabs _anco 
evaluates all the contributions to f|13|) for both scalar one- and two-particle operators. 

The program ANCO can be run in two modes. The interactive mode is intended for testing 
the program and for performing short calculations involving a small number of configurations. 
Normally the modules rabs_anco and rabs_recoupling will be interfaced to GRASP92 or 
to some other program with compatible data structure to perform multiconfiguration or 
configuration-interaction relativistic calculations. The program can handle any combination 
of open subshells with j < 9/2, but subshells with j > 9/2 are only allowed if they contain 
not more than two electrons. 

ANCO is written in Fortran 90/95 and is designed as an addition to the RATIP package \1'A\ . 
The new Fortran 90/95 standard enables us to define new derived data types which will enable 
us to incorporate this module in our continuing development of large-scale computations for 
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open-shell atoms and ions. The full definition of the various derived structures can be found 
in the module header of rabs_anco. Here we need only those types which concern the program 
output. 

All information about pure spin-angular coefficients of scalar one- and two-particle opera- 
tors (t rs (ab) and v^ s (abcd) coefficients) of the Hamiltonian matrix or of some part of it is 
summarized in the derived data anco_pair_list. 

type (anco_csf _pair) , dimension( : ) , pointer :: anco_pair_list 

which is defined by 

type : : anco_csf _pair 

integer : : r, s 

integer :: no_t_coeff, no_v_coeff 

type (anco_t_coef f ) , dimension( : ) , pointer :: t_coeff 
type (anco_v_coef f ) , dimension( : ) , pointer :: v_coeff 
end type anco_csf _pair 

The integers r and s respectively index the bra- and the ket- configuration state functions 
(CSF) for the current matrix element. The variable no_t_coef f is the number of pure spin- 
angular coefficients of one-particle operators that can be constructed for the pair r , s, and the 
variable no_v_coeff gives corresponding data for two-particle matrix elements. The array 
anco_v_coef f contains pure spin-angular coefficients for two-particle scalar operators. It is 
defined by 

type : : anco_v_coeff 

integer : : k 

type(nkappa) : : a, b, c, d 

real (kind=dp) : : v 
end type anco_v_coeff 

where k is the rank k of the effective interaction strength, a, b, c, d point to the relevant 
subshells riiliji, njljjj, niili>ji>, rij/lj/jj/, and the pure spin-angular coefficient itself is given 
in v. The array anco_t_coef f is defined in the same way. 

Memory for the array anco_pair_list is allocated dynamically using the 
allocate ( anco_pair_list (1 :number_of _pair_list_max) ) 

instruction, and can be deallocated subsequently. 

3.2 Interactive calculations 

A typical interactive dialog for calculating spin-angular coefficients is shown in Fig. 1. 
Enter a file name for the anco.sum file: 
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After this prompt, the user should insert the output file name to which the main output 
data will be written. This must be followed by the name of the input file listing CSF in the 
GRASP92 format. The next question 

Generate only not trivial angular coefficients which include (at least one) 
open shells ? 

The response y or Y will cause the program to calculate coefficients for peel shells only; the 
reponse n or N the program calculate will yield data for all shells (open and closed). 

The question 

Generate one-electron angular coefficients for scalar interactions ? 

This requires answer n or N if one-electron coefficients are not needed. If the user responds 
y or Y, then the prompt 

Generate GRASP92-like T coefficients for scalar interactions ?. 

appears. The response y or Y causes GRASP92-like T rs (ab) coefficients to be generated, 
whereas the alternative n or N yields t® s (ab) coefficients. 

A similar dialog follows for two-electron angular coefficients. A number of examples illustrate 
the usage of ANCO in this mode in section 4 below. 

The prompt 

Enter a file name for the anco.vnu file: 

permits the user to specify where the spin-angular coefficients should be stored. 

The module ANCO needs one input file . csl containing the CSF list output by GRASP92 
which is gnerated by the program GENCSL. 

The program creates two output files. The file . sum contains the summary of the problem. 
The other file .vnu contains the angular momentum coefficients and their characteristics. Its 
format for T rs (ab) and V^ s (abcd) coefficients is the same as that of file genmcp.dbg generated 
by the module genmcp of GRASP92 0. The format for t° rs (ab) and v* s (abcd) is very similar 
but without the sorting process which is used in GRASP92. 

3.3 Distribution and installation of the program 

As a new component of the RATIP package |13| similar to the module RCFP ANCO 
will be distributed as a tar archive file of the directory ratip_anco. On a UNIX (or com- 
patible) workstation, the command tar -xvf ratip_anco . tar reconstructs the file struc- 
ture. The directory ratip_anco then contains the Fortran 90/95 modules rabs_anco.f and 
rabs_recoupling . f , the program xanco . f (the main program for interactive work) as well as 
the makefile make-anco. It also includes a number of examples in the subdirectory test-anco 
and a short Read. me which explains further details about the installation. Since the same 
file structure is preserved in both cases, the combination of ANCO with RATIP is simply 
achieved by running the command cp -r ratip_anco/. ratip/. inside the RATIP root 
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ANCO: Calculation of angular coefficients for symmetry— adapted CSF functions 
from the GRASP92 structure program (Fortran 90 version) 
(C) Copyright by G Gediminas and others, Kassel (2000). 

Enter a file name for the anco.sum file: 
test . sum 

Enter the name of the configuration symmetry list file: 
argon- sd. inp 
Loading configuration symmetry list file . . . 

There are 16 relativistic subshells; 

there are 761 relativistic CSFs ; 

. . . load complete . 
Generate only non-trivial angular coefficients which include 
(at least one) open shells ? 

y 

Generate one-electron angular coefficients for scalar interactions ? 

y 

Generate GRASP92-like T coefficients for scalar interactions ? 

y 

Generate two-electron angular coefficients for scalar interactions ? 

y 

Generate GRASP92-like V~k coefficients for scalar interactions ? 

y 

Enter a file name for the anco.vnu file: 
test . vnu 

Figure 1: The typical interactive dialog for calculation spin-angular coefficients. 
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directory; then make -f make-anco will generate the executable xanco, together with the 
other two components xcesd99 |14j and xreos99 13 of the RATIP package. The name of 
the (Fortran 90/95) compiler and special compiler flags can be overwritten in the header of 
the makefile as necessary. Although ANCO uses six other modules which are part already 
of RATIP, no further adaptation of the program is needed. At present, the ANCO program 
has been installed and tested under the Linux and AIX operating systems but, owing to the 
compliance of the Fortran 90/95 standard, no difficulties should arise on any other platform. 

The subdirectory test-anco lists a number of examples which demonstrate the usage of the 
program. 

4 Timing and verification of ANCO 

Tests and timing studies using the Dirac-Coulomb Hamiltonian only were performed for the 
3s 2 3p 6 1 S state of Ar I with the common closed shells ls 2 2s 2 2p 6 for different values of final 
orbital momentum J. The wave function expansions used were: 

1. 3SD: Single and double excitations from 3s 2 3p 6 to the active set {3s, 3p, 3d}contains 14 
configuration state functions (CSF) for J = and 34 CSF (the maximum) for J = 2. 

2. 3SDT: Single, double and triple excitations from 3s 2 3p 6 to the active set {3s, 3p, 3d}. The 
maximum number of CSF is 145 for J = 2. 

3. 4SD: Single and double excitations from 3s 2 3p 6 to the active set 

{3s, 3p, 3d, 4s, 4p, Ad, 4/}. The maximum number of CSF is 465 for J = 2. 

4. 4SDT: Single, double and triple excitations from 3s 2 3p 6 to the active set 
{3s, 3p, 3d, As,Ap,4d, 4/}. 

5. 5SD; Single and double excitations from 3s 2 3p 6 to the active set 
{3s, 3p, 3d, 4s, Ap, 4d, 4/, 5s, hp, 5d, 5/, 5g}. 

We first considered simple cases with a small number of CSF (3SD, 3SDT, 4SD with J = 0, 
1, 2, 3, 4, 5, 6, 7, 8, 9). Although ANCO generates the full set of "pure" coefficients for both 
one- and two-particle operators, the calculation runs from 1.4 - 2.3 times faster than an 
equivalent calculation with GRASP92 because of the lower computational overheads. Tabic |2] 
demonstrates similar enhanced performance for the much larger 4SDT and 5SD eamples, 
showing the improvement expected for large-scale calculations. 

From the results presented in the Table[2]we conclude that the new program is not much faster 
in simple cases, but does better in more complicated cases. The fact that ANCO calculates 
approximately twice the number of angular coefficients as GRASP92 reduces the effective cost 
per coefficient by a further factor of two. 

Although the program is completely new, we have verified that the results presented agree 
completely with those obtained from GRASP92. We have also verified that the Breit inter- 
action is treated correctly, although no data are presented here. 
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Table 2: Timing Comparison for GRASP92 and ANCD codes. Times are given in hours, minutes, 
seconds 





ASF 


Number of 


Running time of 


Speed 


expan. 


CSF 


Tf s (ab) or 
t k rs (ab) 


V r k s (abcd) 


V 


k s (abed) 


GRASP92 


ANCO 


-up 


4SDT 


(J= 


0) 


2 149 


3 


606 


756 023 


1 


530 086 


00:08:11 


00:03:01 


2.7 


4SDT 


(J= 


1) 


5 786 


14 


017 


4 070 156 


8 


188 130 


00:59:01 


00:15:19 


3.9 


4SDT 


(J= 


2) 


8 016 


21 


356 


7 018 885 


14 


077 044 


01:42:47 


00:26:28 


3.9 


4SDT 


(J= 


3) 


8 378 


21 


342 


7 634 136 


15 


290 955 


01:53:55 


00:30:37 


3.7 


4SDT 


(J= 


4) 


7 284 


15 


971 


6 111 074 


12 


260 139 


01:33:17 


00:23:01 


4.1 


4SDT 


(J= 


5) 


5 349 


9 


435 


3 810 165 


7 


656 054 


00:50:27 


00:14:18 


3.6 


4SDT 


(J= 


6) 


3 370 


4 


556 


1 836 602 


3 


706 544 


00:21:52 


00:06:40 


3.3 


4SDT 


(J= 


7) 


1 788 


1 


789 


693 761 


1 


412 443 


00:07:26 


00:02:29 


3.0 


5SD 


(J= 


0) 


468 




621 


75 192 




150 455 


00:00:32 


00:00:17 


1.9 


5SD 


(J= 


1) 


1 134 


2 


324 


395 450 




792 560 


00:03:10 


00:01:29 


2.1 


5SD 


(J= 


2) 


1 609 


3 


704 


697 651 


1 


395 839 


00:06:27 


00:02:44 


2.4 


5SD 


(J= 


3) 


1 584 


3 


441 


721 907 


1 


444 095 


00:06:43 


00:02:59 


2.3 


5SD 


(J= 


4) 


1 361 


2 


500 


558 223 


1 


117 681 


00:05:15 


00:02:17 


2.3 


5SD 


(J= 


5) 


920 


1 


361 


314 909 




632 306 


00:02:30 


00:01:22 


1.8 


5SD 


(J= 


6) 


559 




644 


141 328 




284 102 


00:01:02 


00:00:36 


1.7 


5SD 


(J= 


7) 


259 




226 


44 137 




89 398 


00:00:15 


00:00:12 


1.3 
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5 Examples 



To illustrate the use of ANCO in its interactive mode, we studied Ar II. The program genmcp 
of GRASP92 jHj is used first to generate the CSF file argon-sd. inp before running the 
program xanco. 

The first example does a GRASP92-style calculation. After checking all triangular condi- 
tions for A^, OMiomfe (a/97<5) (see expression (88) in |12]). it multiplies each pure-two particle 
coefficient by the factor X^, oulomb (aP^y5) and prints all non zero V^ s {abcd) coefficients. With 
corrsponding T rs (ab) coefficients, this generates a total of 11279 coefficients. The second 
example, with the same input, calculates a total of 433911 pure non-trivial spin-angular 
coefficients at one go, as is more convenient for large scale calculations. However, only 11279 
coefficients from this set are required in the first example. 

The Test Run Output displays the . sum files and the first 10 lines of . vnu files for both 
examples. 
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TEST RUN OUTPUT 

Example 1 



»type example 1 . sum 
ANCO run at 15:50:39 on Feb 14 2000. 

There are 17 electrons in the cloud 
in 761 relativistic CSFs 
based on 16 relativistic subshells. 
Total number of pair is: 289941 

Generate only not trivial angular coefficients 
there are 716 Grasp92-like T coefficients 
there are 112078 Grasp92-like Vk coefficients 
the total number of coefficients is 112794 



>> type examplel.vnu (first 10 lines of examplel.vnu file) 

ANCO 

761 16 4 



V" 


[( 


2)]_[ 


1, 


1] 


( 3p 


, 3p ; 


3p , 


3p ) = 


-1 


200000000000E-01 


V" 


[( 


1)]_[ 


1, 


1] 


( 3p 


, 3d-; 


3d-, 


3p ) = 


6 


666666666667E-02 


V" 


[( 


3)]_[ 


1, 


1] 


( 3p 


, 3d-; 


3d-, 


3p ) = 


-2 


571428571429E-01 


V" 


[( 


1)]_[ 


2, 


1] 


( 3p 


, 3d ; 


3d-, 


3p ) = 


-1 


632993161855E-01 


V" 


[( 


2)]_[ 


3, 


1] 


( 3p 


, 4s ; 


3p , 


3d-) = 


2 


529822128135E-01 


V" 


[( 


1)]_[ 


3, 


1] 


( 3p 


, 4s ; 


3d-, 


3p ) = 


-2 


108185106779E-01 


T" 


[ 


6, 


1] ( 


4d-, 


3d-) 


1 


000000000000E+00 


V" 


[( 


0)]_[ 


6, 


1] 


( Is 


, 4d-; 


Is , 


3d-) = 


2 


000000000000E+00 



Example 2 

>>type example2.sum 

ANCO run at 16:35:02 on Feb 14 2000. 



There are 17 electrons in the cloud 
in 761 relativistic CSFs 
based on 16 relativistic subshells. 
Total number of pair is: 289941 

Generate only not trivial angular coefficients 
there are 716 pure one-particle angular coefficients 
there are 433195 pure two-particle angular coefficients 
the total number of coefficients is 433911 

>> type example2.vnu (first 10 lines of example2.vnu file) 

ANCO 
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761 16 4 



pure two-particle [ 


( 1)]_[ 


1, 


1] 


( 3p , 


3p ; 


3p , 


3p ) = 


5 


000000000000E- 


-02 


pure two-particle [ 


( 2)]_[ 


1, 


1] 


( 3p , 


3p ; 


3p , 


3p ) = 


-1 


500000000000E- 


-01 


pure two-particle [ 


( 3)]_[ 


1, 


1] 


( 3p , 


3p ; 


3p , 


3p ) = 


5 


000000000000E- 


-02 


pure two-particle [ 


( 1)]_[ 


1, 


1] 


( 3p , 


3d-; 


3p , 


3d-) = 


3 


000000000000E- 


-01 


pure two-particle [ 


( 3)]_[ 


1, 


1] 


( 3p , 


3d-; 


3p , 


3d-) = 


-2 


000000000000E- 


-01 


pure two-particle [ 


( 0)]_[ 


1, 


1] 


( 3p , 


3d-; 


3d-, 


3p ) = 


-2 


500000000000E- 


-01 


pure two-particle [ 


( 1)]_[ 


1, 


1] 


( 3p , 


3d-; 


3d-, 


3p ) = 


-2 


500000000000E- 


-01 


pure two-particle [ 


( 2)]_[ 


1, 


1] 


( 3p , 


3d-; 


3d-, 


3p ) = 


-1 


500000000000E- 


-01 
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